Image transfer through a chaotic channel by intensity correlations 
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The three-wave mixing processes in a second- order nonlinear medium can be used for imaging 
protocols, in which an object field is injected into the nonlinear medium together with a reference 
field and an image field is generated. When the reference field is chaotic, the image field is also chaotic 
and does not carry any information about the object. We show that a clear image of the object be 
extracted from the chaotic image field by measuring the spatial intensity correlations between this 
field and one Fourier component of the reference. We experimentally verify this imaging protocol 
in the case of frequency downconversion. 
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Since the very beginning of its history, second-order (x^) nonlinear optical processes, also called three- wave mixing 
(TWM), have been used to implement image processing. Initially the process involved was frequency upconversion, 
namely second-harmonic generation 0, Q> S S @- Though contem por arily proposed 0, image-processing in 
downconversion took some more time to be realized 0, 0, |g, El DU LLJ- H3 • The imaging capability of the nonlinear 



\ interactions are connected with the spatial properties of the process, which in turn reflect the presence of phase- 

and frequency-matching conditions among the interacting light fields^]. In the past years, we have investigated, both 
theoretically and experimentally, the common features between image processing in TWM and the Gabor holography 
E3, EE EE O , that were first recognized by Firester [3j,y,|5|. As in an holographic process, TWM imaging protocols 
require the presence of a reference field, a role that can be played by any of the interacting fields. The generated image 
thus depends on the simultaneous presence of both the object and the reference fields: this opens the possibility of 
performing conditional measurements exploiting the correlations between the reference and the image fields. Imaging 
Qh, protocols based on correlations have been used in the quantum domain, in which the generation of couples of entangled 
+!> ■ photons by spontaneous downconversion and their detection by coincidence techniques allow to transfer spatial 
information from one of the twin photons to the other (l9l |20L l2lH22| . In these protocols, often called ghost imaging, 

a ' 
cr 



the object is inserted on either signal or idler beyond the crystal. Similar results have been obtained by using a 
classical source of correlated single-photon pairs j2J, 13 ■ In the many photon regime, it has been theoretically shown 
that the same kind of image transfer can be implemented with both quantum entangled [2{| and classically correlated 
light prl I27I l2gj . One may alternatively place the object to be imaged on the beam pumping the spontaneous 
k> , parametric downconversion process |29| . Also in this case, images have been actually detected by mapping suitable 
j_j ■ photon coincidences between the single-photon pairs in the generated beams [3(3, 001 • 

In this paper we exploit the correlation properties intrinsic to \ processes to recover the image of an object 
transmitted through a chaotic channel. In fact, if we use a chaotic reference field for the TWM interaction, we 
obtain the generation of a chaotic image, that does not contain any recognizable information on the original object. 
Nevertheless, we will demonstrate that by evaluating the intensity correlations between one of the Fourier components 
of the reference field and the chaotic generated image, we can recover a clean image of the object. As shown in [33], this 
imaging protocol can be considered a way to simulate many-photon quantum-imaging experiments by using classical 
fields. 

The paper is organized as follows: in Section [H] we present the analytical description of the imaging process 
implemented by a x^ interaction in a nonlinear crystal and and experimentally verify the theoretical predictions for 
a 2-D object; in Section IlIII we show, both analytically and experimentally, how to use the intensity correlations to 
recover the image of the object in the case of a chaotic seed. 
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II. IMAGE TRANSFER IN PARAMETRIC DOWNCONVERSION 



In this Section we present a full 3-D theory of three wave mixing and apply it to the case of a seed field (at ui\) 
interacting with an intense pump field (at oj 3 > u>i). The interaction produces the amplification of seed field and 
the simultaneous generation of a field at u>2 = u> 3 — u>\- This interaction can be used for transferring an amplitude 
modulation from the pump beam to the generated beam. We demonstrated that this imaging process has holographic 
properties by showing that, when the amplitude modulation is set by a 3D object, the wavefronts of the generated 
field reconstruct a 3D image of the object |2J|. In the experimental paragraph of the Section we verified the theoretical 
expectations as to location and sizes of the images. 



A. Theory 

According to the geometry of Fig. ^ we describe the interaction inside a nonlinear crystal of the following three 
amplitude-modulated plane-waves, propagating along kj 



Ei (r,i) 
E 2 (r,i) 
E 3 (r,i) 



x / 2i]ohuji 



X / 2l] hlU2 



"2 



ai (r) exp (— iki • r) + c.c. 



a2 (r) exp (— ik.2 ■ r) + c.c. 



n 3 



[y a 3y (r) + z a 3z (r)] exp (-ik 3 ■ r) + c.c. , 



(1) 



where rij are the refraction indexes, t]q is the vacuum impedance and the optical frequencies, ujj, satisfy energy 
matching (uj 3 = u>2 + In the "slowly varying envelope" approximation, the system describing the second-order 
interaction inside the crystal is 

k x • Vai = i [g+a 3y (r) + g-a 3z (r)] a* 2 (r) exp [-iAk • r] 

k 2 • Va 2 = i [g+a 3y (r) + g^a 3z (r)] a\ (r) exp [-iAk • r] 

k 3 • \7a 3y = ig+ai (r) a 2 (r) exp [iAk ■ r] 

k 3 • Va 3y = ig-ai (r) a 2 (r) exp [iAk ■ r] , (2) 

where g± are coupling coefficients |36j and we have defined the phase mismatch vector as Ak = k 3 — k 2 ki . 

In order to analytically solve J5J), one of the interacting fields must be taken as non evolving during the interaction: 

we are interested in the case of non-evolving transverse field a 3 (r) = a 3 (0). system © becomes 



ki • Vai = ig a 3 (0) a* 2 (r) exp [—iAk • r] 
k-2 ■ Va2 = ig a 3 (0) a* (r) exp [—iAk • r] , 

where g is a new effective interaction coefficient. The solution of Q is given by 



(3) 



ai (r) 



a 2 (r) 



ai (0) 



cosh | Q~y ■ r 



i—— sinh ( Q • r 



(Ak • ki) Q 



exp 



(Ak-ka )Q V 2 / 



O2(0) 



Ak 



, Ak \ .Ak . / Ak 
cosh Q— t smh Q— -r 



exp 




(4) 
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where 




Q=\ ^W^-An . (5) 



^Ak • kj) (Ak ■ k 2 ) 

When the wavevectors satisfy the phase-matching condition (Ak = 0), solutions (0J become: 

n tn , , ( 9\aa (0)1 c \, a 3 (0) . , / g\a 3 (0)| c N 
oi (r) = ai (0) cosh — - b • r + a 2 (0) — - smh — — ^ b • r 



b-k 12 / |«3 (0)| I b ki 
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, x . . / n x gg (0) . , / g|a 3 (0)| C \ , , n , , / ffl«3 (0)} /; \ , R , 

a 2 (r) = » 0l (0) | -— | smh b ■ r + a 2 (0) cosh b • r (6) 

|a3(0J| \bk 12 / \bk 12 / 

where b is along the bisector of ip, the angle ki-to-k 2 , i.e. b = (ki + k 2 )/2 . Since cos?/> can be written as 

cos-0 = sin/3i sin/3 2 + cos/3i cos/32 cos (^i — i? 2 ) , (7) 
we find that b = cos(^/2)(= b • ki j2 ) to be used in ©. Thus (b • r)/(b • ki :2 ) = /(#,/3)r, where 

„ „ s sin/?i + sin/3 2 + cos/3i sini^i + cos/?2 sini? 2 + cos/?i cos??i -(- cos/?2 cos??2 ,„ s 

® = (8) 

is a pure geometrical factor. In the following we will consider the particular case of a 2 (0) = and ai(0) ^ in the 
weak-conversion approximation (<? (03 (0)| /(#, /3) r <C 1), for which the solution © can be written as 

ai (r) ~ ai (0) 

a 2 (r) = iff Mi9)roI (O)oa(O) , (9) 

that is the amplification of field ai is almost negligible and the generated field ai is linearly dependent on both the 
pump field 03 and the incident field a±. 

Note that @ gives the values of the evolved fields at any position r inside the crystal, given a starting point r = 
at the entrance face. Different field components starting at different initial points evolve independently and solution 
© can be thus generalized for any transverse field distribution incident on the crystal. In particular, if field ai(0) 
is a plane wave and the non-evolving pump field 03(0) is amplitude modulated, © shows that the modulation is 
transferred to the generated field a^. 



B. Image formation 



We will now explicitly calculate the "image" field generated on given an "object" field a 3 in a particular 
propagation scheme that will be useful for experimental purposes. 

According to Fig. [21 we consider a modulation of the pump field af(xo,yo) on the object plane (xo,yo)- We 
insert a converging lens (focal length /) on the plane (xl, jjl) located at a distance d from the object plane and at a 



distance dp from the crystal entrance face, plane ( x f , i/f) 
in the presence of an infinite ideal lens is given by |35| 



The field distribution at the entrance face of the crystal 
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By evaluating the second integral and letting do 

fa 



2/ (2/ — 2/ system) we simplify the expression as 



a 3 (t f ) 



■ exp 



.k 3 , 2 



y%) 



2irid 

dx dy a$ (x ,yo) exp 
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(11) 
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where d = 2/ — d F is the distance between the crystal and the image plane of the lens. We now let the field impinging 
on each point of the crystal entrance face interact according to © by setting 0,3 (0) — > a 3 (r F ). According to Q, the 
spatial amplitude of the generated field E 2 becomes 



E 2 (rout) oc a 2 (r out - r F ) exp [-ik 2 • (v out - r F )\ 

oc igrf(d, (3) a* (0) a 3 (r F ) exp [-ik 2 • (r out - r F ) 



(12) 



being z out — zp = L the crystal depth and r = |r otti — vp\. In the limit of non-evolving pump, we can take 
substitute in (|T2"|> 03 (rp) with 03 (r „ t ). Now we let field E 2 freely propagate from plane (x out ,y ut, z out ) to plane 
{ x 2i V2i Z2) along its direction. By neglecting refraction at the crystal interfaces, the field will travel a distance 
s 2 = (z 2 — z out )/ (cos $ 2 cos f3 2 ) (see Fig. 0). By inserting l(TT|) into itT^ and applying free propagation, we get 



E 2 (v 2 ) cx igrf(d,P) k2 h - ^(0)exp 
2ms 2 Lmd 



k 2 (xl+yf\ 



2s 2 



exp [— ik 2 (s 2 + cos j3 2 cos"d 2 L)] exp [ik 2 (sin f3 2 xp + cos/3 2 siri'& 2 yp) 



a% (xo, yo) exp 



fe 3 / -d , 2 2 v 

— — ( x o + yo) 



2d f 

x exp \-ik 2 (s\n.f} 2 x out + cos j3 2 siwd 2 y out )] 



exp 



x exp 



i — (x 2 x out + y 2 y ut) 
s 2 



exp 



. fc 3 + vlut) 



2d 



i—r (xoutxo + youtyo) 

d 



dx ou t dy Q ut • 



dx dy 



(13) 



By exchanging the integrations and neglecting the quadratic term inside the last integral we get 
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£ 2 (r 2 ) oc igrf(#,(3) h h ! , a?(0)exp 
lms 2 Zmd 
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(14) 



where we have introduced the phase factor $ (r 2 ). By writing all the coefficients in a single complex function C we 
can summarize our result as 



E 2 (v 2 ) = Ca* 1 (S))a < i 



k 9 d 



(x 2 - x 2 ) 



k 9 d 



(2/2 - y2. 



,k 3 s 2 ' k 3 s 2 

where we have defined 

x. 2 = s 2 sin (3 2 ; y 2 — s 2 cos (3 2 sin § 2 . 
If we locate the plane (x 2l y 2 , z 2 ) so that s 2 ~ k 2 /k 3 d |38|. we can simplify the expression as 

E 2 (r 2 ) = Cal(0)a^ (x 2 - x 2 ,y 2 - y 2 ) . 



(15) 

(16) 
(17) 



The result in (|17f) shows that with this choice of the propagation configuration, the field distribution on a plane 
satisfying the holographic relation among the distances is equal to the field distribution on the object plane, provided 
a double inversion of the coordinates and a shift in position which depends on the propagation direction Ei. Note 
that the insertion of the imaging lens on the plane (i/,, y^) was necessary to obtain a real image |3q. 
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C. Experiment 



We experimentally verified the imaging properties of the interaction by using the setup depicted in Fig. [3] 
The nonlinear crystal was a type I /3-BaB 2 04 crystal (BBO, cut angle 32°, 10 mm x 10 mm x 4 mm, Fujian Castech 
Crystals). The object field was produced by locating on the plane (xo,Vo) a copper sheet with three holes (~ 256 
/im diameter, see inset of Fig. 0) as the mask producing object O. The imaging lens located in the plane (xl,Ul) at 
a distance do = 60 cm from the object plane had a focal length / = 300 mm so as to realize a 2/ — 2/ system. The 
fields Ei and E3 entering the crystal were obtained from the fundamental (Ai = 1064 nm) and the second harmonics 
(A3 = 532 nm) outputs of a Nd:YAG laser (10 Hz repetition rate, 7-ns pulse duration, Spectra-Physics)thus A 2 = 
1064 nm. The BBO crystal was located at a distance <1f = 20 cm from the lens, so that the image plane of the 
lens resulted to be at a distance d — 40 cm beyond the crystal. A system made of a polarizing beam splitter plus a 
half-wavelength plate was used to obtain a ordinarily polarized field Ei. For the present experiment, the diffuser D 
in Fig.[3]was removed from the setup. The sensor of a CCD camera (Dalsa CA-D1-256T, 16 /im x 16 /im pixel area, 
12 bits resolution, operated in progressive scan mode) was located in the detection plane [X2, yi-,z-i)- The distance of 
the CCD from BBO was chosen to be s 2 = 20 cm, so as to satisfy s 2 = d k 2 /k 3 . In Fig. 0] we show the image of the 
holes as detected by the CCD camera. Note that in agreement with H17|) the transverse dimensions of the image were 
equal to those of the object. 



In this Section we consider the same interaction seeded by a chaotic field Ei, which amounts to realize the same 
image transfer process of above through a chaotic channel. We will demonstrate that in this "chaotic" situation the 
imaging system gives no straightforward information on the object. Nevertheless, the intensity fluctuations intrinsic 
to the chaotic light can be profitably used to implement a protocol for correlated imaging that allows the recovery of 
the object image. 



We now suppose that the seed field Ei in Q can be written as an incoherent superposition N plane waves having 
random complex amplitudes, ai, n , and wave vectors, ki iM , with random directions but equal amplitudes, fci, n = 2n/Xi 



According to solution JSJ, field Ei does not evolve inside the crystal so that its spatial part can be written as 
E\ {T ou t) — E\ {tf). Due to the chaotic nature of field Ei, we expect that also its Fourier transform is chaotically 
distributed. In fact, if we insert a lens of focal length / on the path of field Ei at a distance d from the exit plane of 
the nonlinear crystal, so as to have the Fourier plane on {x%, yx,Zi) we get 



III. CHAOTIC IMAGE-TRANSFER IN PARAMETRIC DOWNCONVERSION 



A. Theory 




(18) 




(19) 
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By substituting the definition of Ei 



JY 



Ei(ri) oc ^ 

ai ' n 2irif 6XP 



. fei (x 2 +y 2 ) ( _ t 



2/ 



/ exp (-iki >n • r out ) exp i-i (xix onf + yiy ou t) 



dx ut dyout 



N 



53 C^ai^exp 



2/ I /, 



' k\ ~ fei 

x <5 ( -j(xi - /sin^i^), - /cos/3i )n sini?i,„) 



(20) 



which means that if the amplitudes a\. n have random values, the intensity distribution on the plane has the form of 
a speckle field: 



N 

h (n) oc \Ei (ri)| 2 = 53 |C"| 2 \ a i,n\ 2 S (xi - fsmpi tn ,yi - f cos /3 i)n cos i?i in J , 



(21) 



Inside the nonlinear medium, each of the spatial Fourier components of the seed field that is phase matched with the 
pump field generates an independent contribution to E2 according to @ . The overall field at the crystal output face 
is given by 



N 



E 2 (rout) oc ^2igr n f(d n ,{3 n )a* l n (0)a 3 (r F )exp[-i k 2j „ • (r out -r F )] , 



(22) 



where each wavevectors k.2 irl are assumed to satisfies the phase-matching conditions k.3 = ki .„ + k2, n - As obvious, 
this assumption cannot be satisfied for all wavevectors \s-i >n , thus setting an angular limitation to the effectiveness of 
the interaction. The result (fTK|l can thus be used to evaluate each of the terms in l(2*2*|l to find the N contributions to 
field E 2 



E2 ( r 2,n) = C n a\ n a 3 



&2 d ,_ k 2 d 

X2.n - X 2 ) , — (V2.n ~ V2 , 



k 3 s 2 . 



k 3 S 2 ,n 



where we have defined 



x 2 ,n = s 2 ,„sin/3 2 ,„ ; y 2 , n = s 2 , n cos f3 2 , n sinz9 2i . 



(23) 



(24) 



If the effective angular spread of the wavevectors k 2 „ is not too broad, we can assume that s 2j „ ~ dfc 2 //c3 for all 
components and that the corresponding N images form on the same plane. Thus from (|24|l we can write 



E 2 (r 2 ) = 5Z Cn a l,n fl 3 ( X 2,n ~ X 2 ,y 2 ,n ~ V2) ■ 



(25) 



n=l 



Note that also the the field on plane (x 2 , y 2 , z<£) results to be the sum of N holographic images of the pump field E3, 
reconstructed in N different transverse locations. If we now take into account the random nature of the amplitudes 
di,n, we obtain that the sum in l|25|l is as incoherent as that in l|18[l . We can thus write the detected intensity of field 
E2 in the plane (x 2 , y 2 , z 2 ) as 



JY 



I 2 (r) oc \E 2 (r 2 )| 2 = 53 \ C n\ 2 |ai : „| 2 \a 3 (x 2 ,„ - x 2 ,y 2 ,„ - y 2 )\ 



(26) 



n=l 



This result means that no information about the object field is any more detectable in a direct way. Neverthe- 
less, we can use the correlation properties between fields Ei and E 2 introduced by the nonlinear interaction in 
the crystal to recover the image of the object. To do this, we select a single position (xi,yi) in the plane of 
the Fourier transform of Ei, that, according to l|2ip. corresponds to choosing a single plane wave having ampli- 
tude a±j. We then calculate the correlations of the temporal intensity fluctuations between field Ei and field E 2 we get 
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(A/iA/ 2 ) = {hh) - (h)(h) 

N 

= (\ a lj\ 2 ^\Cn\ 2 \ai, n \ 2 \a% (X2,n ~ X2,V2,n ~ ll2)\ ) 
n=l 
N 

- {\ a lj\ 2 )(^\Cn\ 2 \ai,n\ 2 \a*i (X2,n ~ X2,V2,n ~ V2)\ ) 
n=l 

JV 

= ^2(\ a l,3\ 2 \ C n\ 2 \ai,n\ 2 \a£ (X2,n ~ X2,V2,n - V2)\ > 
n=l 
N 

~ ^2(\ a l,j\ 2 )(\ C n\ 2 \ai,n\ 2 \a° (X 2 ,n ~ X 2 ,y 2 . n ~ 2/2 ) | } 

n=l 

JV 

= ^2 l^™' 2 1°? ( i2 ^« ~ x 2,V2,n ~ 2/2 ) | 2 ((Kjl 2 |ai : „| 2 } - (|ai,j| 2 )(K,n| 2 )) 
n=l 
N 

= ^2 \C„\ 2 \af {X 2 ,n - X2,y2,n ~ IJ2) f & 2 (|ai )T ,| 2 J 5j. n 

71 = 1 

= IQlVfW,/) \a$(x 2J -x 2 ,y 2 , j -y 2 )\ 2 , (27) 

where a 2 (x) =< x 2 > — < x > 2 is the variance. In deriving (|27H we used the fact that the coefficients C n do 
not depend on the amplitudes ai,„ (see l|14[) and QI5[0 and that the pump field amplitude 03 at each sample of the 
statistical ensemble. The result in ()27|l shows that the correlation function between the image on field E 2 and a single 
point on the Fourier plane of field E\ can reconstruct the image of the object encoded on field E3. 

B. Experiment 

For the experimental verification of the results of the previous section, we used the same setup in Fig. [31 modified 
by introducing the light diffuser D on the seed beam E\. The diffuser, a ground-glass wheel, was moved from shot 
to shot of the laser in order to obtain the temporal statistics needed to evaluate the correlation function H27JI . A 
portion of the chaotic field Ei emerging from the diffuser was selected with an iris (PH, in Fig.[3J| of ~ 8 mm diameter 
and then filtered in polarization with a polarizing beam splitter, PBS, and a half-wave plate. Lens L2 (/ = 15 cm) 
provides the Fourier transform of Ei on the plane (x\,y\, z\). The detection planes of (x\,y±, Z\) and (x 2 ,y 2 , z 2 ) were 
made to coincide on the sensor of the same CCD camera so that each signal occupies half sensor. 
First of all, we checked the chaotic nature of both the seed field E\ and the generated field E 2 both in space and in time. 
A characteristic of chaotic light is to have an intensity obeying thermal distribution, P(I) = exp(// < I >)/ < I >. 
We thus measured the spatial probability distributions, P r (h) and P r (/ 2 ), of the intensity recorded by the different 
CCD pixels for a single shot, relative to the Fourier transform of the seed field E% and to the generated field £'2- 
Figure [S] shows the results for the spatial distributions along with the intensity maps used to evaluate them (insets). 
From these results we can see that both the spatial distributions can be fitted by a thermal distribution and that 
the intensity map relative to the generated field E 2 has no memory of the field modulation (the same three holes as 
before) imposed on the pump. 

We also measured the temporal (over many repetitions of the laser pulse) probability distributions, Pt(Ii) and 
Pt(h), of the intensity recorded by choosing a single pixel of the CCD in the map of E\ and E 2 and recording the 
intensity values at each laser shot. Also in this case the intensity probability distributions are well fitted by thermal 
distributions (see Fig.^J. The temporal traces of the intensities of the selected pixels are shown in the insets of the 
figure. 

Once established the correspondence of our experimental setup with the requirements of the theory, we evaluated the 
correlation function l|27|l over 1000 shots by taking the whole map I 2 (x 2 ,y 2 ) and by selecting the value of a single 
pixel in the intensity map of I\ . 

In Fig. [3 we show the resulting reconstructed image (map of G (hj, I 2 {x 2 ,y 2 ) j): the similarity in the quality of this 
image compared with that obtained without diffuser (see Fig.0J is really impressive, in particular if the reconstructed 
image is compared with any of the single-shot intensity maps I 2 {x 2 , y 2 ), see for instance the inset of Fig- EI 



9 



IV. CONCLUSIONS 

In conclusion, we have demonstrated that the spatial intensity correlation properties of the downconversion process 
can be used to recover a selected image from a chaotic ensemble of holographic images. Note that theoretical results 
similar to those in Section Til Al and Al would be found for any choice of the object and reference fields among the 
three interacting ones. The image recovered by G (Ix j, ^2(^2, IJ2)) fulfils the holographic properties of the difference- 
frequency generated hologram that would be obtained by using the single plane- wave Ei.j as the seed/reference field. 
We expect that the method also works in the case of an unseeded process, in which no reference field enters the crystal. 
In this case any twin beam in the parametric fluorescence cone would play the role of our reference and image fields 
and our intensity correlation protocol should provide an a posteriori selection of a single holographic image. Note 
that the plane on which this recovered image forms would depend on the choice of the plane-wave component of the 
twin party used as the reference. 
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FIG. 1: Interaction inside the crystal. The shaded plane contains the optical axis and the normal to the crystal entrance face 
(2— axis). 
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FIG. 3: Experimental setup: HS, harmonic separator; D, diffusing plate (removable); Mi_s, mirrors. Lens Li images O into O' 
through a 2/ — 2/ system. Lens L2 realizes the Fourier transform of field E\ on the CCD sensor. Inset: image taken with an 
optical microscope of the copper sheet containing the three holes used to produce the object-field modulation. The diameter 
of the holes was ~ 200 /im. 
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FIG. 4: Single-shot holographic image detected in the plane (£2,2/2,22) 
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FIG. 5: Probability distributions, P T (Ii,2) of the intensity recorded by the different CCD pixels for a single shot, relative to 
the Fourier transform of the seed field E\ and to the generated field i?2- Insets: single-shot intensity maps of h and h- 
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/ (arbitrary units) 



FIG. 6: Probability distributions, Pt{I\,2) of the intensity of a selected CCD pixel recorded for many successive laser shots. 
Insets: temporal traces of the intensities. 



